o 
o 
o 

Switching of Geometric Phase in Degenerate 
(M : Systems 

■ _ A. Yahalom'^'* and R. Englman ^''^ 

^ ' " Faculty of Engineering, 

Tel-Aviv University, Tel-Aviv 69978, Israel 
Research Institute, 
College of Judea and Samaria, Ariel 44284, Israel 
^ • Department of Physics an Applied Mathematics, 

Soreq NRC,Yavne 81810, Israel 
^ [ e-mail: asya@post.tau.ac.il, englman@vms.huji.ac.il 

O ; February 1, 2008 



C/3 



, Abstract 

I The geometric and open path phases of a four-state system subject 

to time varying cyclic potentials are computed from the Schrodinger 
] equation. Fast oscillations are found in the non-adiabatic case. For 

' parameter values such that the system possesses degenerate levels, the 

, geometric phase becomes anomalous, undergoing a sign switch. A 

' physical system to which the results apply is a molecular dimer with 

two interacting electrons. Additionally, the sudden switching of the ge- 
ometric phase promises to be an efficient control in two-qubit quantum 
computing. 

C3 ■ PACS: 03.65. Ge; 03.67.Lx; 76.60.-k 

o ■ 
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X : 1 Introduction 



Following the discovery of the topological (or geometric) phase acquired by 
a system that evolves adiabatically around a closed path in a parameter 
space |l|]-[^, many instances of the phenomenon have been found either 
theoretically or in experiments A widespread quantum mechanical 
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application of the theory is when an electron in a molecule slowly evolves 
subject to a field due to (say) the host ions, such that the field varies peri- 
odically in time. In typical cases the electronic state changes sign after a full 
period [p. It is entirely obvious that for a pair of electrons, which are inde- 
pendent except that both are subject to fields with the same period (though 
the fields may have otherwise different characters) the total electronic wave 
function (being a product of two individual wave functions) will return to its 
original value after a period. Similarly, when the field on the second electron 
is a constant or has a period which is ^ or | (and so on) times the period for 
the first field then, after one complete revolution, the total wave function 
will again acquire a change of sign. It is of interest to consider next what 
happens to the geometrical phase if one introduces an interaction between 
the two electrons. One would expect that the geometrical phase would not 
be altered by the interaction. This is indeed the case for the model stud- 
ied in this work, with one notable exception: When the two-electron states 
become degenerate by some coincidence of the parameters in the Hamil- 
tonian, then under suitable conditions the expected sign change becomes 
reversed. This phenomenon is studied in this work. Other aspects of states 
that have degeneracies on or off the trajectories have been noted earlier in 
p. The topological phase change after a full revolution has been inter- 
preted in several previous publication as a surface integral in the parameter 
space H, 1^. We find the acquired phase by a different method: namely, by 
numerical integration of the time dependent Schrodinger equation, which is 
the appropriate method for circumstances that the surface integral is not 
readily available. Moreover, this method, which was previously used for a 
smg le electron 0, 0, |l2|, gives the entire open path phase, rather then the 
geometrical phase alone. Theoretical expressions and interpretations for the 



open path phase were given in O, 14]. The last section contains algebraic 



and numerically less convenient calculations. 



2 The physical model 



Brief descriptions of the model have been given previously 16]. In 
essence, we consider a dimer consisting of two vibrating, planar and pos- 
sibly dissimilar molecules lying upon each other. A simple instance is two 
triatomic molecules ^3 - B3. In addition to several tightly bound electrons, 
each molecule contains one electron that is in an orbital doublet state, de- 
generate in the symmetric (equilateral triangle) configuration of its host. 
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The molecules are sufficiently remote for the orbital doublets to be local- 
ized on the host molecule and also for overlap (and exchange) effects to 
be negligible. There remain vibronic interactions between the electronic 
and the in-plane vibrational motions, as well as an inter-electronic coupling. 
The former (Jahn- Teller) interaction |17| is conveniently treated by the for- 



malism of Longuet-Higgins et al [y^], which represents the two electronic 
states by e^*^™ , where 6^ an angular variable and m = 1,2 for the two 
molecules. In this work we shall modify this formalism by using as the ba- 
sis the direct product of two real representations, namely cos cos etc. 
The non-totally symmetric vibrational mode coordinates of the molecules 
are denoted by qm cos (pm and Qm sin . In terms of the electronic and 

vibrational coordinates, as given above, and the vibronic coupling strengths 
Gm, the interaction may be simply expressed as 

Hm = Gm cos(26'm - (pm) m = 1, 2 (1) 

^Prom elementary physical considerations, the vibronic constants Gm are 
(at least approximately) proportional to the mode amplitudes Qm- The 
interaction between electrons on the two molecular entities is familiar from 
energy transfer and other subjects, and is expressed by terms that contain 
jointly excitation and de-excitation operators on two molecules ]l9|. For 
simplicity we assume for the interaction term the large U [t/U 0)," anti- 
ferromagnetic Heisenberg" limit of the well-known Hubbard Hamiltonian 
[pO| , having the product form 



Hi2 = 2G cos(26ii - 0i) cos(26i2 - ^2) (2) 

This is symmetric and even in the relative angular variables and is invari- 
ant under a full rotation of either mode angular variable <j)m- Analogously 
to Hubbard's U which gives the exchange splitting between like and unlike 
spins on the same site, G represents the electrostatic splitting between like 
and unlike vibronic states on the two molecules. G depends foremost on 
the distance between the molecules and only weakly on the molecular dis- 
placement coordinates. In previous works [^, the cases Gm larger and 
smaller than G were studied separately. (To exhibit formally the weakness of 
overlap effects, one would introduce, both in the wave functions and in the 
couplings, two more factors that depend on two further coordinates, (say) ri 
and r2, such that the overlaps between the ri and r2 factors are negligible.) 
The total Hamiltonian H is the sum of equations (|^) and (^). 

H = Hi + H2 + H12 (3) 
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We next represent the vibrations of the molecules as rotational motions 



in the {q, <f>) -planes. (The potential surfaces in Figure 3.3 of |]17[ admit 
of such motion.) Classically, this is equivalent to taking qi = constant, 
q2 = constant, (pi = toit, (/>2 = u)2t, where t is time and uji and U2 are the 
angular frequencies of the (small) molecular displacements. The equation 
that forms the basis of this work is the time dependent Schrodinger equation, 
with fi = \, 

'a = if* (4) 

The two-electronic product states are labeled as follows: 

tpi = N cos 9i cos 02 

11)2 = N sin 6i cos 62 

ips = N cos 61 sin 62 

= iV sin 6*1 sin 6*2 (5) 

where the normalization constant, N = -. In this basis the Hamiltonian 
finds the following representation: 

2 

/Gi cos + G2 COS 1^2 C?i sin 01 + G sin 01 cos (1)2 G2 sin 02 + G cos 0i sin 02 G sin 0i sin 02 

-\-G cos 01 COS 02 

Gi sin 01 + C? sin 0i cos 02 ~Gi cos 0i -f G2 cos 02 G sin 0i sin 02 G2 sin 02 — G cos 0i sin 02 

— G cos 01 cos 02 

G2 sin 02 + G cos 01 sin 02 G sin 0i sin 02 Gi cos 0i — G2 cos 02 Gi sin 0i — G sin 0i cos 02 

— G cos 01 cos 02 

VG sin 01 sin 02 G2 sin 02 — G cos 0i sin 02 Gi sin 0i — G sin 0i cos 02 — Gi cos 0i — G2 cos 02 

+ G cos 01 COS 02 

(6) 

The eigen energies of this matrix, denoted by K^^{r = 1,..,4) and the 
corresponding eigen functions (the so-called adiabatic eigen states) are: 



4"" = i(Gi + G2 + G) ^ i;^"" = - cos(0i - ^) cos(02 - 

4"" =\{-Gi + G2-G) ^V2^^ = isin(0i-^)cos(02-$) 

11 z z 

i(Gi -G2-G) ^ i^i"" = - cos(0i - ^) sin(02 - ^) 

TT Z Z 



'^3 — 2 
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i(-Gi-G2 + G) ^vr = -sin(0i-^)sin(02-^) (7) 

TT Z Z 
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3 Time evolution 



Recalling that (pi = uit and (j)2 = ^2^-, we substitute the time dependent 
Hamiltonian matrix, equation (^) , into equation (Q) to obtain the solu- 
tion ^ (which is now a four column-vector) by forward integration, using a 
" Mathematica" numerical algorithm. With the choice of wi = 1 and u:2 a 
simple multiple of wi, the periodicity is 27r and the solutions exhibited in the 
following figures stretch over a single period of the Hamiltonian, t = Q — 2-k. 
Beyond this, the solutions can be obtained by symmetry. The initial value of 
^' at t = was taken as (1, 0, 0, 0). This choice implies no loss of generality, 
since any other solution can be derived from this through a change of the 
parameters, shift of the time scale and superposition of solutions. The ana- 
lytic treatment of the Schrodinger equation including the interaction, given 
in the last section , involves the solution of a quartic equation which is not 
decomposable into two quadratic equations. This stands in contrast to the 
one-electron two- state problem, whose Schrodinger equation involves the 
solution of a quadratic [21, 22, 12 1. Thus the interaction introduces in the 



problem a genuinely new element. 
(A) Non-adiabatic effects. 

We first investigate the approach to the adiabatic limit, which is the 
condition under which the results of Q were derived. To achieve this limit, 
at least one of the coupling coefficients Gi, G2 and G has to be large. 
(Basically, their ratio to omega serves as the measure of adiabaticity |11]. 
We define it by A = |^|) The approach to adiabaticity is exemplified in 
Figure 1, whose two curves have identical parameters except that all coupling 
strengths are reduced by a factor 10 in the thinly drawn curve , thus not 
yet reaching the adiabatic limit. 

In either curve there is a sign change in the wave function at full revolu- 
tion, but the lack of adiabaticity shows up in two ways: First, in the state 
amplitudes it induces oscillations, whose periods depend on the coupling 
strength. These oscillations are the ionic-field analogues of the Rabi oscilla- 
tions (that are due to excitations by an electro- magnetic field). The former 
have been analyzed in detail in an earlier paper on a two-state model jl2U . 
Secondly, the state amplitude does not reach -1, but only approximately 
-.75. This reduction has already been quantitatively explained for the two 
state model in |10| and will not be further discussed here. In the fully adia- 
batic limit, the solution can be read off from the analytic expressions given 
in equation (0) . Having started (by virtue of the initial conditions imposed) 
with ipi^ , we see that the thick curve shown in Fig. 1 indeed approximates 



5 



Re Amplitude 
1 

. 5 



-0 . 5 
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Figure 1: Non-adiabaticity effects in tlie real part of the initially excited 
component, as a function of time. The frequencies on the two dimers are 
iOi = 1 and UJ2 = 2. The values of the coupling parameters are as follows. 
Thick line: Gi = —80, G2 = —160, G = 40 (near adiabatic limit). Thin 
line: Gi = —8, G2 = — 16, G = 4 (non-adiabatic case) 

. The dynamic phase has been subtracted in this and in the following figures. 



well to the expression cos(t/2) cos(i), which is the amplitude of the first 
component in ipf'^-, according to equation (^) and equation (|^ . Further 
solutions for the near-adiabatic limit are shown graphically in Figures 2 and 
3. 

(B) Enters degeneracy. 

Further computations in the adiabatic limit (large G's) with other pa- 
rameters likewise confirm the expressions in equation (|^) , as predicted by 
the adiabatic theorem. As is well known, this theorem is based on the non- 
crossing of energy levels ||2^ . In the absence of inter-electronic coupling, the 
two electrons move independently and, therefore, even if there was a degen- 
eracy in the system through the combination of the levels, no effect would 
be felt. For a coupled two-electron system, the degeneracy is expected to 
produce new effects. 

This was indeed found. Thus, when in equation (^) the inter-electron 
coupling was so adjusted that 

G = -Gi (8) 

then the first two levels in equation (0) become degenerate. The accom- 
panying curves (a)-(c) in Figure 4 show the way the amplitude (of the first 
component) changes as this degeneracy is traversed through subsequent se- 
quence of changes in G. The results are similar when the changes are made 
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Figure 2: Two-electron state amplitudes in a dimer. The thick line shows 
the time dependent amplitude of the first (initially excited component), the 
thin line that of the second component in Equation (5). Frequencies: loi = 1, 
UJ2 = 4, Gi = —40, G2 = —80, G = 16 (near adiabatic limit) 

ReAmplitude 




time 



Figure 3: Two-electron state amplitude in a dimer, with both molecules 
subject to a periodic force. After a full revolution the two electronic states 
each change their sign, leaving the total state invariant. Frequencies: lui = 1, 
= 1, d = -100, G2 = -200, G = 40 (near adiabatic limit). Thick line: 
First, initially excited component. Medium thick line: Second and third 
components. Thin line: Fourth component. 
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Figure 4: Degeneracy and near-degeneracy effects in the first, initially 
excited component. The following values are common to the three curves: 
uji = 1,L02 = = —40,^2 = —80. The inter-electronic coupling G is 

varied as follows: Thin- line: G = 36 and 44 (non-degenerate). Medium 
thick line:G = 39.6 and 40.4 (near degenerate case). Thick line: G = 40 
(exact degeneracy). 

in an electron- vibration coupling strength Gi. The curves are computed in 
the near-adiabatic limit. 

One observes that with degenerate parameters, the amplitude keeps its 
sign at a half period and achieves at full revolution a switched sign with 
respect to that in the non-degenerate situation. The switch occurs abruptly 
as the coupling strength(s) is (are) varied. We now turn to a discussion of 
this switch, regarding its dependence on adiabaticity and its stability, an 
application and the algebraic treatment of the results. 

4 Discussion 

The preceding figure (fig. 4) was plotted in the near adiabatic limit when 
A = is large(= 40). In the extreme adiabatic limit when ^4 = 00 only the 
symmetric ( thin ) line (ending at -|-1) and the antisymmetric (thick) line 
(ending at —1) are possible, the latter occurring at exact degeneracy and 
the former for any values of the parameters that do not cause exact degen- 
eracy. When the adiabaticity parameter A is finite, though large, for large 
deviation from degeneracy the computed curve approaches the symmetric 
curve and for small values of the deviation , the antisymmetric curve. It is 
possible to describe this effect quantitatively, by finding a relation between 
the adiabaticity A and that value of the deviation from precise degeneracy 
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at which the real part of the ampUtude at full revolution (at time = 2tt) is 
zero (i.e., just half way between the values far from and exactly at degener- 
acy) . Attaching to this value of the deviation from degeneracy the subscript 
0, we find from our computation for the present model 



(We recall that at G + Gi = one has a strict degeneracy.) For a numerical 
illustration of this result, when the adiabaticity j4 = 40 (this being the abso- 
lute value of the ratio of the coupling strength Gi to the angular frequency 
of motion uJi, as e.g. when Gi = — 40,u;i = 1), then a deviation from de- 
generacy less than 0.54 will cause in the phase a switch whose magnitude 
is at least ^ .In fig. 4, the thin line is for a deviation 4 that is larger than 
0.54 and the phase switch is close to zero , the medium thick line is for a 
deviation of .4, which is below the 0.54 limit and the phase switch is larger 
than ^.Though equation (^) is derived from numerical fitting (of over one 
and a half orders of magnitude in A and to an average accuracy of about 
five percent) for a specific model, we expect that a qualitatively similar re- 
lation would hold in other cases where the degeneracy causes a switch in the 
phase. A deviation from (or removal of) the degeneracy can also achieved 
by adding off-diagonal terms 5Hi2 = 6H21 in the (1,2) and (2,1) positions 
to he Hamiltonian matrix of equation(6).When these perturbation terms are 
constants their effect on the switching is given very closely by equation (^) 
, when the deviation on the left hand side stands for SH12, or the added 
matrix element. This represents a test for the stability of the switch against 
a small deviation from the special form of the interaction assumed in this 
work. (It would be tempting to combine the two findings into a single for- 
mula, in which the left hand side of equation (^) is simply the splitting of 
the two near degenerate levels, given by 



but our numerical results indicate that the above formula holds without the 
factor 4.) 

A further stability test was also made: When a constant was added to 
each of the terms in top-right to bottom-left diagonal, the computed curves 
were not changed in a perceptible way, even when the constant reached .2 
of the (root mean square) value of the term. 



|(G + Gi)o| =0.54^ 



(9) 



splitting = {{G + Gi)l + 4{6Hi2f)2 



(10) 
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A vibrating dimer, carrying one electron on either molecule, was used as 
a model system to introduce the formalism. We now turn to an alternative 
application. 

The four state formalism in this paper can be naturally interpreted as a 
two-qubit system. "Qubits" are two-state building-blocks for quantum com- 
puting. (An introductory account is found in ||2j].)The two-qubit situation 
is the minimal one to describe logical operations between control and target 
[ p5| , 26], and it has also been shown to be sufficient for that purpose |27]. 
For an identification of logical operations in a two-qubit scheme, one can 



refer to |25]. The correspondence between the time dependent description 



(or the Schrodinger equation @), on which the present work is based) and 
quantum computing has been developed in |29|. One of the several possible 
implementation is the use of trapped ions in doublet quantum states jSOtj . 
as a control-no gate. This may be done as follows. We have seen that as 
the first and second components in equation (0) are brought to degener- 
acy through the choice G — > —Gi, then (upon full revolution) there is a 
sign-switch in the first component in equation (0) . Likewise, in the second 
(co-degenerate) component. There is, however, no switch in the third or in 
the fourth components (which are energetically distinct). Suitable combi- 
nations (e.g., sum and difference) of the first two components can thus be 
used for representation of "Yes" and "No" operations, with the choice of the 
coupling parameters acting as control. One important result of this paper 
is the abrupt switching of the geometric phase through a slight change in 
the coupling strength, as shown in Figure 4, and this can be conveniently 
used for manipulation on the target, or of the outcome, by external control. 
The abruptness of the switching is of some additional interest, in that the 
change in the total phase of the target, upon a slight change of the param- 
eter, is more due to the geometric phase than to the change in the dynamic 
phase. This means that the dynamic phase need not necessarily be sub- 
tracted from the total phase. Abrupt switching of the Pan char atnam phase 
for two polarized light beams was described and observed earlier 

Further extension of our results might involve consideration of several 
(more than two) qubit states. This may be applicable to quantum computa- 
tion, in which an important issue is how to achieve reduction of decoherence 
of phases by random causes. It has thus been proposed that a possible way 
is the use of superposition states composed of several qubits 24 1. 
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5 Algebraic treatment 



The solutions of the Schrodinger equation (^) and equation (P) , obtained in 
the text numericahy, can be exhibited through a sequence of(three) matrix 
multiplications operating on a column (four) - vector which contains 
the initial conditions. The vector is written in the representation of equation 
(6) as 

/ e 



""^*Xi(0) \ 



(11) 



*X4(0) / 

and contains initial values needed to yield the initial values of the solution 
and the eigen frequencies Kr (of the Rabi-type) of a four matrix. In the 
extreme adiabatic limit and excluding degeneracy, the eigen frequencies 
can be given explicitly and are identical with the eigenvalues shown in 
equation (0) . Non-adiabatic corrections to the (non-degenerate) eigenvalues 
are of the order of (^)^ . The full solution is 



(12) 



The matrix M generates a complex representation when applied to the 
real basis shown in equation dSf) . Its form is 



( 1 
1 
1 

V 1 



-1\ 



-1 



(13) 



The next matrix T is diagonal and contains time exponents with the coef- 
ficient ^, which is the source of the sign change upon a full revolution. 



/ g-5(wi+<^2)i 



J" 



\ 









62 










62 







|(-a;i+a;2)t 







62 








|(a;i+a;2)t 



(14) 



Lastly, A diagonalizes the following matrix with the eigenvalues Ky. (r 
1,..4) appearing above. 
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H 



St 





/ 


—Ul — UJ2 


Gi 


G2 


^ \ 


1 




Gi 


LOl — UJ2 


G 


G2 


2 




G2 


G 


—UJl + UJ2 


Gi 




I 


G 


G2 


Gi 


UJl + UJ2 / 



(15) 



substitute into the Schrodinger equation 



To prove equation equation ([ 

(4) and transpose the time derivative of the factor in equation (14) from 
the left to the right hand side. In the adiabatic hmit, when the frequencies 
are vanishingly small compared to the coupling strengths, A simplifies to 



/ 1 1 1 1 \ 

1-1 1 -1 

1 1 -1-1 

V 1 -1 -1 I J 



(16) 



When some of the adiabatic eigenvalues approach degeneracy, e.g.. 



.AD 



AD 



(17) 



the situation changes dramatically. The corresponding non-adiabatic eigen- 
values differ from the adiabatic ones by a term that is linear in a frequency. 
Specifically, to a first order in the w's. 



.AD 



2 



(s = l,2) 



(18) 



This half-frequency value in the k -exponents shown in equation (|ll| ) adds 
on a further factor of —1 upon a full revolution, thus accounting for the 
sign switch found in the degenerate situation. Also (as the numerical so- 
lutions show), two components in the four- vector become vanishingly small 
throughout the time-range. With the initial conditions chosen for our solu- 
tion, the vanishing components are the second and the fourth. If we remove 
the corresponding rows and columns from the Hamiltonian shown in equa- 
tion @ , the resulting two-by-two matrix becomes under the degeneracy 
condition Gi = — G, 



H 



1 [ (G2 + Gcos 0i) cos ( 

2 I {G2 + Gcos(l)i)sm( 

scalar) 



{G2 + Gcos <pi) sin(/)2 
— {G2 + G cos (pi ) cos (p2 



1 



-Gcos 
2 ^ 



a 



(19) 
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Since {G2 + Gcos(/>i) is a common factor for all elements in the matrix, it 
does not enter the component amplitudes and these contain only the angle 
<j)2- This represents a second algebraic justification for the phase switch 
upon degeneracy. 

The relation in equation (p|) represents a three dimensional surface of 



degeneracy in the four dimensional displacement-manifold [15, 16 1. Actu- 
ally, the surface is (as expected) only two dimensional, for the following 
reason: The degeneracy involves states 1 and 2. An off-diagonal coupling 
between these states would lift the degeneracy, even if this coupling would 
be so weak as not to affect the dynamic behavior. By our ignoring this cou- 
pling, we effectively put it equal to zero, which represents a further relation 
between some of the coordinates, say, between qi and (72 -In summary, we 
have investigated both the geometric phase (i.e., at full revolution) and the 
interim open path phases for two interacting electrons, each in a doublet 
state, when they move in cyclically varying potentials. With the choice of a 
product inter electron potential, the eigen energies are constants along the 
motion. The dynamic phase (being the time integral of an eigen energy) 
can thus be conveniently subtracted off from the total phase. With more 
general interactions, one would need a more complicated procedure for sub- 
traction. We have identified non-adiabatic effects (namely, high frequency 
oscillations and the incomplete return to the initial state after a full revolu- 
tion). We then elucidated a switch of the sign of the wave- function as the 
system develops a degeneracy in the two-electron energy levels. The switch 
was shown to be quite sudden as the degeneracy situation was approached 
through varying the parameters in the Hamiltonian. 
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